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We present results from a fully relativistic Monte Carlo simulation of diffusive shock acceleration (DSA) in 
unmodified (i.e., test-particle) shocks. The computer code uses a single algorithmic sequence to smoothly span 
the range from nonrelativistic speeds to fully relativistic shocks of arbitrary obliquity, providing a powerful 
consistency check. We show the dependence of the particle spectrum, which can differ strongly from the 
canonical f(p) oc p~ 423 result, on the obliquity and on the strength and "fineness" of scattering for a range of 
shock Lorentz factors. The Monte Carlo results are also shown to be consistent with a simple relation for the 
spectral index given by Keshet and Waxman| 1 1 in parallel shocks when the diffusion is "fine" and isotropic. 

1. Introduction 

Most work on diffusive shock acceleration (DSA) in relativistic shocks has been restricted to the ultra-relativ- 
istic regime with magnetic fields assumed parallel to the shock normal. However, trans-relativistic shocks are 
certain to be important in some sources and, in general, relativistic shocks will have highly oblique magnetic 
fields. Here we consider test-particle diffusive shock acceleration in shocks of arbitrary obliquity 1 and with 
arbitrary Lorentz factors, 70. We consider only particle acceleration in plane, unmodified (i.e., test-particle) 
shocks where effects on the shock structure from superthermal particles are ignored. In contrast to nonrelativ- 
istic shocks, the details of particle scattering strongly influence the superthermal particle populations produced 
in trans-relativistic and ultra-relativistic shocks (see 1 1 1|2|| 12 1 for example). Unfortunately, these details are 
not known with any reliability so the results of all current models of DSA depend on the particular scattering 
assumptions made. We use a simple, parameterized model of particle diffusion (described in detail in .6*1 ) 
which we believe adequately describes the main features of particle transport. Until a self-consistent theory of 
wave-particle interactions in relativistic shocks is produced, parameterization will be necessary. 

2. Model and Results 

The acceleration process is modeled with a Monte Carlo simulation where thermal particles are injected far 
upstream and followed as they diffuse through the shock (for full details, see and references therein). 
Scattering is elastic and isotropic in the local plasma frame and the scattering mean free path A in the local 
frame is proportional to the gyroradius r g . We simulate small-angle scattering by allowing the tip of the 
particle's fluid-frame momentum vector p to undergo a random walk on the surface of a sphere. After a small 
time St the momentum undergoes a small change in direction of magnitude 56 within a maximum angle <5# m ax- 

If the time in the local frame required to accumulate deflections of the order of 90° is identified with the 
collision time t c = X/v p , |7| showed that (5# m ax = ^J&St/t c , where St is the time between pitch-angle 
scatterings. We take A proportional to the gyroradius r g = pc/ (eB) (e is the electronic charge and B is the 
local uniform magnetic field in Gaussian units), i.e., A = ?7 m fp r g , where 77 m f p is a measure of the "strength" of 
scattering. The strong scattering limit, ?7 m fp = 1, is called Bohm diffusion and in this limit, cross-field diffusion 
is important. Setting St — r g /N g , where N g is the number of gyro-time segments St, dividing a gyro-period 
T g = 2-KTgjvp, we have S9 max = ^12ir/ (r) m [ p N g ). Large values of N g result in "fine" scattering where it 
is assumed that magnetic field fluctuations exist on all scales so that 89 can be arbitrarily small. The DSA 



1 Oblique shocks are those where the angle between the upstream magnetic field and the shock normal, #bo < i s greater than 0° . Parallel 
shocks are those with 0bo = 0° . Every where in this paper we use the subscript (2) to indicate upstream (downstream) quantities. 
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distribution, f[p), becomes independent of N g when it is large enough such that 59 <C l/7o> otherwise large- 
angle scattering (LAS) effects become important. Thus, the scattering properties of the medium are modeled 
with the two parameters 77 m f p and N g and, along with 70 and 6*30, this yields, in the high Mach number limit, 
four independent parameters characterizing DSA. 
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Figure 1. Test-particle spectra for shock Lorentz factors 70 = 3 (left panel) and 70 = 10 (right panel) (note that p 4 ' 23 f{p) 
is plotted). For 70 = 3, the shock compression ratio, as determined in 1 41, is r = 3.16. For 70 = 10, r = 3.02. The listing 
of line types applies to both panels. The light-weight solid curves labeled LAS show the effects of large-angle scattering, 
all other curves have 89 <C 1/70- All spectra are calculated at the shock in the shock frame. 

In Fig.^ we show examples where these four parameters are varied. We plot p 423 f(p) for a trans-relativistic 
70 = 3 and a more fully relativistic 70 = 10. In all cases except the two marked LAS, N g is large enough so 
that 86 <C l/7o- In the two LAS cases, the effects of large-angle scattering are clearly seen. 

In the limits 59 <C I/70 and 70 3> 1, DSA produces a power law f(p) oc p~ 4 23 , independent of 9bo or 
f?mfp (see 0). The heavy solid curve in the right panel of Fig. [2 shows that for 9bo = 0° and rj m { p = 1, the 
canonical power law is obtained. However, even with 59 <C I/70, f(p) can be much steeper than this for larger 
values of 9bo and 77 m f p . We showed in [6 1 that effects of 9bo and r) m f p diminish with increasing 70 but persist 
noticeably to at least 70 = 30. 

The curves labeled LAS in Fig.[2both have 9bo = 0° and r/ m f p = 1, but here N g is small enough so particles 
make large deflections in individual scattering events. Large-angle scattering will be important whenever 
TJmfpNg is small enough so that 59 max <jt. I/70 and, in principle, could occur even in the shocks expected 
when a 7-ray burst fireball expands ultra-relativistically into the interstellar environment. 

Keshet and Waxman[10| have presented the following analytic expression for the power-law spectral index 

ctkw 

<tkw = (3/3 - WoPl + Pi) I (fa - 2 ) , (1) 

where (3q (fa) is the speed of the upstream (downstream) fluid normalized to the speed of light. This expression, 
which assumes that the shock is parallel and that 59 <C I/70, depends only on the shock Lorentz factor and 
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Figure 2. Comparison of Monte Carlo spectra (solid histograms) with the analytic result given in 1 10 1 (dotted lines). 
p 4 ' 23 f(p) is plotted and the normalization is arbitrary. 



the compression ratio. The compression ratio must be determined independently. In Fig.|2]we compare Monte 
Carlo spectra to eq. Q for several 70. The Monte Carlo results use the solution to the shock jump conditions 
(including the equation of state) calculated in [4 1 to determine r and this value is used in eq. Q (i.e., 02 — flo/ r) 
and listed in Fig. |2] The match between the two models is quite good in this trans-relativistic range and will 
improve outside of this range. 



3. Discussion and Conclusions 

In nonrelativistic shocks, the test-particle power law from DSA depends only on the shock compression; it is 
independent of the details of diffusion, on the strength of scattering, and on the obliquity (i.e., independent of 
N g , T] m fp, and 6bo)- In relativistic shocks, particularly trans-relativistic ones, test-particle spectra can depend 
strongly on these parameters. We have shown examples with various N g , r) m f p , and 9&o and, in general, the 
power law steepens if rj m i p > 1 and/or 9bo > 0°. The parameter N g determines the fineness of scattering and 
if N g is small enough so that 69 max 5^ l/7o> large departures from a power law can occur and the spectrum 
can become harder than p~ 4 23 even for ultra-relativistic shocks. In addition, as have shown that the simple 
analytic result given in eq. Q 1 10 1 matches our Monte Carlo results for parallel shocks and when 89 <C 1 /7o- 

The parameterization of diffusion used in the Monte Carlo technique is useful for investigating the behavior 
of DSA in relativistic shocks. It makes fewer approximations than analytic models and can cover a far larger 
dynamic range than any current particle-in-cell (PIC) simulation. Nevertheless, the actual diffusion depends 
on the self-consistent generation of magnetic turbulence by the accelerated particles and is inherently complex 
and has not been described adequately even for nonrelativistic shocks (see, for example, II II ). In fact, the 
correct form for the magnetic power spectrum in self-generated, relativistic turbulence will not be known until 
PIC simulations can perform this calculation. Unfortunately, the treatment of these problems using PIC codes 
is beyond current computing capabilities for three basic reasons: (i) To correctly describe cross-field diffusion 
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and other physics of the viscous subshock, PIC simulations must be done fully in 3-D |8||9|. If 1- or 2-D 
simulations are used, cross-field diffusion, an essential element in DSA, is unrealistically suppressed; (ii) If 
electrons are to be understood, electrons and protons must be modeled simultaneously (except for pair plasma 
shocks), and the simulation must treat widely disparate inertial scales, greatly adding to the run time; and (iii) 
In order for nonlinear effects to become apparent or field generation to occur on large scales, simulations must 
be run long enough in a large enough box with enough particles for a significant population of extremely ener- 
getic particles to be produced. The combination of these three requirements means that, while PIC simulations 
can add to our knowledge of critical aspects of the injection problem and the start of magnetic field generation 
process (e.g., 1 13 1), they will not be able to model the injection and nonlinear acceleration of electrons and ions 
to energies relevant to ultra-high energy cosmic rays or GRBs until computers more powerful than exist today 
become available. 

While we have only discussed test-particle acceleration, nonlinear effects, where the accelerated particles mod- 
ify the shock structure, may be important in relativistic shocks (see [5 1), and in fact, most models of 7-ray bursts 
assume that relativistic shocks efficiently accelerate electrons. If this is the case, test-particle spectra may not 
be good approximations. This is certainly true for DSA in trans-relativistic shocks, where the harder intrinsic 
power law will naturally lead to a nonlinear feedback on the shock structure. 

Nonlinear effects will also influence the distribution of energy between electrons and protons (e.g., |3|) and 
modify the shape of the distribution at low energies, i.e., near the sharp peaks at ~ 70 in Figs. ^ an d|2](i.e., 
1 5 1). These issues are critical for 7-ray burst models. 
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